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^ i' Abstract. It has recently been suggested that observed galaxy rotation curves can 

■^ i be accounted for by general relativity without recourse to dark-matter halos. Good 

fits have been produced to observed galatic rotation curves using this model. We show 
^SJ ' that the implied total mass is infinite, adding to the evidence opposing the hypothesis. 
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1. Introduction 

Understanding tlie nature and origin of the apparent dark matter responsible for the 
bulk of the gravitational mass of galaxies and clusters is one of the most important 
challenges in modern cosmology. Although many candidates have been proposed for the 
required weakly interacting cold dark matter l2|, to date there has been no convincing 
detection of a suitable particle candidate. For this reason there is considerable interest in 
a recent proposal jS| that accurate modeling of galactic dynamics using general relativity 
is consistent with observed fiat rotation curves and luminous mass without the need for 
additional dark matter. In this scenario, the authors assume a cylindrical stationary 
space-time containing pressureless dust and solve the weak-field Einstein equation to 
first order in the gravitational constant G to obtain the associated velocity distribution 
as described below. 

At first glance, this would seem to be a compelling solution to the dark matter 
problem. However, objections have already been raised by several authors, on a various 
grounds. These include the conditions at z = [UEj, general approximation principles 
jHj, and specific analysis of the approximations used [Zj. The original authors have 
responded to some of these criticisms |Sj ■ In this paper, we suggest yet another possible 
problem with this interpretation of the dark-matter problem. We reanalyze the galactic 
rotation curve as in Pj and find that the model, as applied, contains an implied large 
(infinite) mass distribution extending beyond the observed galactic disk, and hence at 
best has simply moved the dark-matter to the outer Galaxy, and more likely, the model 
is ruled out on energetic grounds. 

2. The Cooperstock and Tieu Model 

Our goal is to analyze the implied mass distribution in the Cooperstock & Tieu j2j model. 
In order to deduce this we first summarize their approach. In their analysis of rotation 
curves, the galaxy to be described is treated as a uniformly rotating axisymmetric system 
containing pressureless fluid. The general stationary axially symmetric metric is chosen 
of the form 

ds^ = -e''-'"{udz^ + dr^) - r^e-'"d(j)^ + e"'(crft - Nd(j)f , (1) 

where u, u, w and A^ are functions of cylindrical polar coordinates r, z, and the metric 
coefficient u is taken to be unity. In a co-moving frame the the fluid four- velocity is just 

U' = Si, . (2) 

Making a local (r, z held fixed) transformation to 

(j) = (j) + uj{r,z)t , (3) 

diagonalizes the metric and, in the weak field limit, leads to simple approximations for 

the local angular velocity u and the tangential velocity V: 

Nee"" Nc ,, ,,, 

'^ = -^ iTP^ — ~ ^r ) V = Lur . (4) 



General Relativistic Galaxy Rotation Curves: Implications for Dark Matter Distributions 

The Einstein field equations to first order in the gravitational constant G are easily 
combined to yield a Poisson-like equation for u. 

V^„ + ^.?^, (5) 



where 



V^W = Wrr + Wzz H ' ■ (6) 

r 
The geodesic equation 

rlTP 

— + Tl,U'U' = (7) 

together with the adopted metric and four-velocity (PJ) and (j21) simply becomes 

Wr = Wz = , (8) 

and hence 

V^^ = , (9) 

within the fluid. The interior equations for p and A^ then become: 

(10) 

(11) 







N„. + iV,, - 


r 






Nj + Nl 


8nGp 






r2 


c2 


Eq. 


m 


can be expressed 
V2$ = 


as 


where 










*./> 





(12) 



(13) 

Hence, flat-space harmonic functions $ become the generators of the axially symmetric 
stationary pressure- free weak fields. Using Q and (fT!?|) . the tangential velocity 
distribution becomes 

V = c- (14) 

r 

Since the field equation for p is non-linear, the galactic rotation is deduced by first 
finding the required generating potential $ and from this the appropriate function N 
for the galaxy is deduced. Then Eq. (fTTjl yields the density distribution. 

Each galaxy analyzed requires its own composing elements to build the generating 
potential. A separation of variables yields a solution for $ in terms of Bessel functions 
Jn{kr) 

<l> = Ce-^^'^Jo{kr) (16) 
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Figure 1. Velocity V (m/s) as a function of radius r {kpc) 

where C is an arbitrary constant. The general solution is expressed as a linear 
superposition 

<l> = ^ae-'="I^IJo(W (17) 

n 

with n ~ 10 chosen appropriately for the desired level of accuracy. From ()17|1 and (Q, 
the tangential velocity is 

y = -c5]fc„C„e-^"l^lji(A;„r) . (18) 

n 

This expression was used in jH] to fit the rotation curve for the Milky Way. In 
this paper we analyze this fit and show that any fit implies a large mass at larger radii. 
Hence, this approach does not solve the dark matter problem, but has merely moved to 
dark matter to large radii. Worse, we now have an excess mass problem! 

3. Interpretation of the Model Fits 

Figure HI reproduces the velocity curve fit for the Milky Way given in [Hj, but extended to 
larger radius. The large oscillations in velocity although alarming, are rendered harmless 
if the density is sufficiently low. Figure |21 shows the extended density profile. It would 
be tempting to assume from this that the density could be cutoff at the limits of the 
observed galaxy. However, we show that for any velocity fit the total mass is infinite, 
in contrast to the original claim of finite integrated mass. In fact, smoothed out, the 
enclosed mass rises linearly with radius. 



4. Integrating the mass 

To obtain the total implied galactic mass, consider a model with a single value of k, 

N = re-'''Ji{kr) . (19) 

Eq (fTT|) written for density is 

c2 fN^ + N^\ 
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(20) 
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Figure 2. Density p (/cgm ^) as a function of radius r {kpc) 




Figure 3. rp as a function of radius r (kpc), z = 



In this case 

p = e-2'=^(Jo(A:r)2 + Ji(A:rf) • (21) 

We note rp —^ (2/7r) exp {—2z} as r ^ oo, illustrated in Figure El The mass element is 

dm = pdrdz rdO . (22) 

So the mass enclosed within cylindrical coordinate r is 

M{r) = 27r f drdzrp (23) 

^ 27r('-') /'rf;2e-^" r . (24) 

In this case M{r) clearly rises linearly. 

Now we look at the general case. Taking our cue from Eq.(j2SI), rp is plotted against 
r for the Milky Way model using two different scales of extended range, shown in Figures 
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Figure 4. rp as a function of radius r (kpc) , z = 



4-10 



3-10" 



-p 



a 2-10 



1-10 



y 



m 



u 



u 



u 



y 



1- 



200 400 600 800 1000 

r 



Figure 5. rp as a function of radius r (kpc), z — 



m El It appears that rp oscillates around a positive average value. If the z dependence 
were independent of r like that in the single k case, this would imply that M{r) tends 
to rise linearly in the limit of large r. The following analysis shows this is the case, even 
though the z dependence is not simple. 

If we can find a, (3 > such that for large enough r 
/ dzrp{r, z) 



rp{r, 0) 



> a 



and 



rp{r, 0) > /3r 



then M{r) > a(3r for large enough r. To find a, for fixed r define 

rp{r, z) 



f{z) 



rp{r, 0) 



(25) 
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From PH). 

fiz) = iJ2^ke-'r + iT.he-'^y^ (28) 

k k 

with coefficients ak,bk determined by the modeL Note /(O) = 1. Define a as the 
minimum over all sets a^ybk, in any model, of the integral, 

a= min f dzf{z). (29) 

a > since otherwise f{z) > and continuity of f{z) imply /(O) =0, a contradiction. 



To find j3 first find the form of the density for z = using (f^IH) , 



p(r,0) 



Y^ (rakJoikr) + CkJi{kr)) + ij2'^^kJi{kr) 

\ k J \ k ) 



lr\ (30) 

(31) 
(32) 



with coefficients determined by the model. So as r ^ cxo 
rp^r \^akJQ{kr)\ ^\Y^bkJ\ikr)\ 

-^ I ^ flfc sin fcr j + I XI ^fc cos kr \ , 

with a re-scaling of the coefficients. So the integral 

/ drrp = y^ dr{ak sin krY + (6^ cos fcr)^ (33) 

J ^ Jo 

rr 

"*" X! / '^'" ('^'^1 '^'=2 sin fcir sin k2r + fo^i ^^ cos kir cos ^2^) 

The first sum consists of a piece linear in r plus an oscillating piece, while the second 
sum is purely oscillating. The oscillations are all bounded, so overall the required /? > 

exists, and M{r-) behaves as predicted. 

5. Conclusion 

Despite the close fit to observed velocity curves over the galatic range, we find that the 
total enclosed mass in the model of Cooperstock and Tieu rises as a linearly increasing 
function of radius. This indicates that the model has not removed the dark matter 
problem, but has merely moved the mass to large radius. Moreover, the implied infinite 
mass, and large velocity fields suggest that this interpretation is probably not reahstic. 
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FG02-95-ER40934. 
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